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Abstract. We present results of simulations of stellar collapse and explosions in spherical symmetry for progenitor stars in the 
8-10 Mq range with an O-Ne-Mg core. The simulations were continued until nearly one second after core bounce and were 
performed with the Prometheus/Vertex code with a variable Eddington factor solver for the neutrino transport, including a 
state-of-the-art treatment of neutrino-matter interactions. Particular effort was made to implement nuclear burning and electron 
capture rates with sufficient accuracy to ensure a smooth continuation, without transients, from the progenitor evolution to 
core collapse. Using two different nuclear equations of state (EoSs), a soft version of the Lattimer & Swesty EoS and the 
significantly stiffer Wolff & Hillebrandt EoS, we found no prompt explosions, but instead delayed explosions, powered by 
neutrino heating and the neutrino-driven baryonic wind which sets in about 200 ms after bounce. The models eject little nickel 
(< O.OlSMo), explode with an energy of ^ 0.1 x 10^* erg, and leave behind neutron stars (NSs) with a baryonic mass near 
1.36 Mq. Different from previous models of such explosions, the ejecta during the first second have a proton-to-baryon ratio 
of Te ^ 0.46, which suggests a chemical composition that is not in conflict with galactic abundances. No low-entropy matter 
with Te 0.5 is ejected. This excludes such explosions as sites of a low-entropy r-process. The low explosion energy and 
nucleosynthetic implications are compatible with the observed properties of the Crab supernova, and the small nickel mass 
supports the possibility that our models explain some subluminous Type II-P supernovae. 
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1. Introduction 

Recent observations of subluminous Type II-P supernovae 
(SNe) like 2005cs, 2003gd, 1999br and 1997D, have re¬ 
newed attention to stars near the lower end of the mass range 
of core-collapse SN progenitors, i.e. to stars with about 8- 
10 Mq, which develop O-Ne-Mg cores. A possible link be¬ 
tween both has been suggested because of the low ^®Ni and 
ejecta masses and low progenitor luminosities (e.g., Chugai 
& Utrobin 2000, Hendry et al. 2005). However, due to many 
uncertainties this connection is far from being clear (e.g., 
Pastorello et al. 2004,2005; Hamuy 2003; Zampieri et al. 2003, 
and refs, therein). Also the Crab Nebula’s progenitor was pro¬ 
posed to be in this mass window (Gott et al. 1970, Arnett 1975, 
Woosley et al. 1980, Hillebrandt 1982). The observed composi¬ 
tion of the Crab remnant (small C and O abundances. He over¬ 
abundance) was interpreted as a strong indication that the Crab 
Nebula comes from a collapsing and exploding progenitor with 
an O-Ne-Mg core (Davidson et al. 1982, Nomoto et al. 1982, 
Nomoto 1983). 

Moreover, these stars were considered as possible sites for a 
low entropy r-process (for example, Hillebrandt 1978, Wheeler 
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et al. 1997, Sumiyoshi et al. 2001, Wanajo et al. 2003) based on 
the assumption that they explode by the prompt bounce-shock 
mechanism, which Hillebrandt et al. (1984) found to work in a 
numerical simulation, taking Nomoto’s O-Ne-Mg core model 
(Nomoto 1984,1987). Such explosions are characterised by the 
direct propagation of the shock out of the core, the formation 
of a mass cut, and the continuous acceleration of the material 
outside of the mass cut to high velocities. They are expected to 
eject relatively large amounts of neutron-rich matter with low 
Te (~0.2) and low entropies (~ 10 A:b per nucleon). However, 
several groups could not confirm the viability of the prompt 
explosion mechanism (Burrows & Lattimer 1985, Baron et 
al. 1987, Mayle & Wilson 1988). Mayle & Wilson (1988) con¬ 
tinued their simulations in the post-bounce phase for a longer 
time and obtained instead a so-called neutrino-driven, delayed 
explosion (Bethe & Wilson 1985) with a low production of 
^®Ni (approximately 0.002 Mq) in agreement with subluminous 
Type II-P SNe as mentioned above, but with a vast overproduc¬ 
tion of neutron-rich material (at least 0.02 Mq of ejecta with 
Ye ^ 0.41). The latter finding is inconsistent with the chemical 
composition of our galaxy, which allows for no more than 10“^ 
Mq of material with Ye < 0.42 being ejected per SN (Hartmann 
et al. 1985). Moreover, the explosion energies of both studies, 
around 2 x 10^'erg in Hillebrandt et al.’s (1984) model and be- 
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tween 0.6 and 1.2 x 10^'erg in Mayle & Wilson’s (1988) simu¬ 
lations, would be inconsistent with the long plateau phase of the 
above mentioned subluminous SNe, if their H-envelope masses 
were $ 8 Mo. 

It was suggested that the reason for the discrepant results 
in the SN simulations of O-Ne-Mg cores (prompt explosions, 
delayed explosions, no explosions) could be explained by the 
different nuclear EoSs used by the groups (Fryer et al. 1999). 
Having in mind that the different approximations in the neu¬ 
trino transport in previous calculations introduced additional 
uncertainties, we revisit this topic with a state-of-the-art neu¬ 
trino transport treatment together with a careful description of 
weak interactions and including relevant nuclear burning re¬ 
actions. We additionally make a comparison of collapse and 
post-bounce calculations with different nuclear EoSs. 

2. Numerical techniques and input physics 

The transport of neutrinos and antineutrinos of all flavors is 
done with the energy-dependent solver for the coupled set of 
moments equations and Boltzmann equation called Vertex. It 
is described in detail in Rampp & Janka (2002). The equa¬ 
tions of hydrodynamics are integrated with the Newtonian 
hnite-volume code Prometheus, which uses a third-order, time- 
explicit Godunov scheme. This code is a direct implementation 
of the Piecewise Parabolic Method (PPM), based on a Riemann 
solver. General relativistic gravity is taken into account ap¬ 
proximately by an “effective relativistic potential” according to 
Marek et al. (2005). Gravitational redshift and time dilation ef¬ 
fects are included in the neutrino transport (see Rampp & Janka 
2002). 

The code is augmented with improved microphysics as de¬ 
scribed in Buras et al. (2005). It includes also the improved 
treatment of electron captures on a large variety of nuclei in nu¬ 
clear statistical equilibrium (NSE), based on shell model Monte 
Carlo calculations, as described by Langanke et al. (2003). 
In addition, electron captures on certain important nuclei in 
the non-NSE regime, in particular ^"Ne and ^"'Mg, are imple¬ 
mented according to Takahara et al. (1989). A simplified treat¬ 
ment of nuclear burning accounts for the main reactions of 
seven symmetric nuclei (He, C, O, Ne, Mg, Si, Ni). Details 
about the implemented microphysics will be described in a 
forthcoming paper (Kitaura et al., in preparation). The nu¬ 
clear burning reactions considered by Hillebrandt et al. (1984) 
('^C-h'^C, '® 0 -i-'® 0 , '^Ch-'® 0 ) are all included, taking into ac¬ 
count different reaction channels. 

To describe matter in NSE, we use two different nuclear 
EoSs in separate simulations, the Wolff & Hillebrandt (W&H) 
EoS (Hillebrandt et al. 1984), which is based on Hartree Fock 
calculations, and the Lattimer & Swesty (L&S) EoS (1991), 
which is a hnite-temperature compressible liquid-drop model 
and has a compressibility modulus of 180 MeV. 

This permits us to compare our models with those of 
Hillebrandt et al. (1984), in which the Wolff and Hillebrandt 
EoS was used and which produced prompt explosions. The 
low-temperature and low-density EoS outside of the NSE 
regime is described by an ideal gas of nuclei and nucleons, 
electrons, positrons, and photons (Janka 1999). The switch be- 



Fig. 1. Mass trajectories for the simulation with the W&H EoS as a 
function of post-bounce time (tpb). Also plotted: shock position (thick 
solid line starting at time zero and rising to the upper right comer), 
gain radius (thin dashed line), and neutrinospheres (ve: thick solid; 
Vei thick dashed; v^, v^, Vr'- thick dash-dotted). In addition, the 

composition interfaces are plotted with different bold, labelled lines: 
the inner boundaries of the O-Ne-Mg layer at ~0.77 Mq, of the C-O 
layer at ~1.26 Mq, and of the He layer at 1.3769 Mq. The two dot¬ 
ted lines represent the mass shells where the mass spacing between 
the plotted trajectories changes. An equidistant spacing of 5 x 10 “^Mq 
was chosen up to 1.35794/q, between that value and 1.37654/q it was 
1.3 X lO^^Afo, and 8 x IO^^Mq outside. 



Fig. 2. Velocity profiles vs. enclosed mass at different times for the 
model with the W&H EoS. Times are normalized to core bounce. 

tween NSE and non-NSE description was made in a density- 
and temperature-dependent manner. The initial model is the 
same as the one used in previous SN calculations of O-Ne-Mg 
cores by Hillebrandt et al. (1984) and Mayle & Wilson (1988). 
It is a 2.2 Mq He core that corresponds to a progenitor with a 
main sequence mass of ~8.8Mo (Nomoto 1984, 1987). Prior to 
collapse it has an O-Ne-Mg core with ~13Mq, surrounded by 
a C-O shell of about O.O 8 M 0 . We take, however, the initial data 
at a time when the central density is ^4x10*" g/cm^ and only 
~ 0.1 solar masses at the center of the core have reached nuclear 
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Fig. 3. Velocity profiles as functions of radius for different post¬ 
bounce times for the simulation with the W&H EoS. The insert shows 
the velocity profile vs. enclosed mass at the end of our simulation. 



tpb [ms] 

Fig. 4. Explosion energies as functions of time for the simulations 
with the W&H and L&S EoSs. The energy is defined as the volume 
integral of the total gas energy (internal plus kinetic plus gravitational) 
in regions where the latter is positive. 

statistical equilibrium (NSE). This is earlier than the starting 
configuration taken by other groups, where the core had al¬ 
ready a central density of ~3xl0" g/cm^ and where around 
0.3 solar masses were in NSE. Our earlier initial model allows 
us to trace the evolution of the core towards collapse. 

We added a helium atmosphere of about 10 around 
the O-Ne-Mg and C-O core, so that we could move the outer 
boundary of our Eulerian grid from the core radius of about 
1100 km to 100 000 km. Eor the He-shell profile we adopted a 
power-law like behaviour of the temperature {T oc r ') from a 



Fig. 5. Profiles of the electron fraction Te and entropy s as functions 
of enclosed mass at the end of our simulation with the W&H EoS. 

10.2 Mq progenitor of A. Heger (private communication), and 
constructed the density profile by assuming hydrostatic equi¬ 
librium, a mass fraction of 100% He, and using the EoS for the 
low-density regime. We employ in our simulations a very fine 
mesh in order to resolve the steep density gradient at the outer 
boundary of the C-O layer, with 1600 nonequidistant zones for 
the hydrodynamics part. The neutrino transport is done with 
235 nonequidistant radial zones. 

3. Results 

Our results show a very smooth transition from the preceding 
progenitor evolution to core-collapse. Since O-Ne-Mg cores 
are gravitationally less bound than more massive stellar pro¬ 
genitors and can release more energy due to nuclear burning, 
a temperature- and density-dependent treatment of all relevant 
nuclear burning reactions had to be included, combined with 
a detailed description of the important electron capture rates. 
Only that ensured that the progenitor evolution continued to¬ 
wards gravitational collapse without numerical transients. We 
could therefore confirm that the neutrinos produced by electron 
captures carry away efficiently the energy that is released by 
nuclear burning (Miyaji et al. 1980, Miyaji & Nomoto 1987, 
Hashimoto & Nomoto 1993). A cruder burning treatment, or 
omission of the improved electron capture rates on nuclei, can 
have the consequence that the core expands instead of collaps¬ 
ing to a NS, as we verified in test calculations. Eig. [TJshows 
how the mass shells in the inner region during the first mil¬ 
liseconds start contracting towards the center. The collapse of 
the core proceeds to higher central densities and when the den¬ 
sity of nuclear matter is reached, the EoS “stiffens”, and the in- 
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ner homologous core bounces. The supersonically falling outer 
layers collide with this central core and a hydrodynamic shock 
is formed as can be seen in Fig. 1, where this discontinuity 
in the fluid flow becomes visible by sharp kinks of the mass 
shell trajectories. This happens after 59 ms of collapse for the 
calculation with the W&H EoS and after 78 ms in the cal¬ 
culation with the softer L&S EoS. The shock formation ra¬ 
dius is at 1.40x10® cm, corresponding to an enclosed mass of 
0.475Mo in the former case, and at 1.15x10® cm with an en¬ 
closed mass of 0.425Mo in the second model. This so-called 
prompt bounce shock produces initially positive velocities in 
the postshock matter. However, the energy of the shock is in¬ 
sufficient to cause a prompt explosion, and the photodisintegra¬ 
tion of nuclei consumes such amounts of energy that the shock 
is quickly damped and that within only ~1.2 ms after shock 
formation the velocities are negative everywhere (see Eig.|2for 
the case of the W&H EoS). Therefore, the prompt shock mech¬ 
anism fails, independent of the employed nuclear EoS. 

The subsequent expansion of the shock is supported by a 
combination of different effects. Initially very high mass accre¬ 
tion rates cause the material to pile up between neutrinosphere 
and shock as it is also observed in the early post-bounce accre¬ 
tion phase of more massive progenitor stars (see, for example, 
Buras et al. 2005). Second, the rapid decrease of the mass ac¬ 
cretion rate contributes to ensure ongoing expansion, because 
even for quasi-stationary conditions the accretion shock adjusts 
to a larger radius for smaller mass accretion rates. Einally, as 
soon as the shock reaches the outer edge of the C-O shell, a very 
steep density decline leads to an outward acceleration of the 
shock. The last two aspects are linked to the specific structure 
of O-Ne-Mg cores and discriminate SN progenitors with such 
cores from more massive stars. However, despite the shock ex¬ 
pansion the material behind the shock has initially still negative 
velocities and is accreted onto the forming NS (Eigs. lllSTl . Note 
that when the matter right behind the shock starts to expand 
with the shock, the gas accreted by the shock is gravitationally 
bound and remains so in passing through the shock. However, 
pdV work excerted from below and to a minor extent energy 
input by neutrino heating can convert the accretion into an ex¬ 
plosion, accelerating a tiny amount of matter (« 10“^ Mq) to 
move outward with the shock. 

While the shock reaches larger radii, the temperature and 
density behind the shock decrease. High-energy electron neu¬ 
trinos and antineutrinos, which stream off from their neutri- 
nospheres (represented by the thick solid, dashed, and dash- 
dotted lines in Eig.[3 begin to deposit energy behind the shock 
mainly by absorption on nucleons. This leads to the formation 
of a “gain radius” (thin dashed line in Eig. ^ which separates 
a layer of neutrino cooling around the neutrinosphere from the 
energy “gain layer” behind the shock (Bethe & Wilson 1985). 

About 80 ms after bounce for the L&S EoS and about 60 ms 
for the W&H EoS, the neutrino heating timescale, defined by 
the total energy in the gain layer divided by the neutrino heating 
rate in that region, gets smaller than the advection timescale. 
The latter is given as the time the accreted matter needs for 
being advected from the shock to the gain radius. Since the in¬ 
creasing shock radius leads to smaller and smaller postshock 
velocities, the duration of the deposition of energy via neutrino 


absorption in the shocked matter increases. The continuous in¬ 
put of energy raises the total energy of the gain layer to a value 
near zero within roughly 100 ms, unbinding the matter in the 
gravitational field of the forming NS. The fluid velocity in the 
layer close to the gain radius therefore starts to become positive 
(see Eigs. CD and the explosion energy shows a rapid increase 
(Fig.D. 

The cooling region then becomes more and more narrow as 
the gain radius retreats towards the neutrinosphere, so that neu¬ 
trinos diffusing out of the contracting protoneutron star begin 
to heat the layers right above the neutrinosphere. Gas is thus 
ablated from the NS surface, and the so-called neutrino-driven 
wind phase sets in at f ^ 200 ms after bounce. Several of the 
mass shells depicted in Eig.^clearly show this process. 

At this time the energy in the expanding postshock matter 
has increased to about 0.1 x 10®' erg, rising further due to the 
power input by the neutrino-driven wind. We extrapolate that 
the final energy of the explosion will be slightly larger than 
0.1 X 10®' erg for the calculation with the W&H EoS and might 
be about 50% higher in case of the L&S EoS. This is roughly 
a factor of 10 lower than the canonical SN value, in contrast to 
the findings in previous explosion models of O-Ne-Mg cores 
(Hillebrandt et al. 1984, Mayle & Wilson 1988). 

At the end of our simulation the mass cut, and therefore 
the baryon mass of the protoneutron star, is around 1.360 Mq 
for the W&S EoS (Eig.© and about 1.363 Mq for the L&S 
EoS. The NS mass will only slightly decrease further because 
of the ongoing mass loss in the neutrino-driven wind. The 
mass of the ejecta lies therefore between 0.014 and 0.017 Mq. 
The ejected gas has an electron fraction, Te, between 0.46 and 
0.53 and entropy values between 10 and 40 per nucleon for 
both EoSs (see Eig.|3for the case with the W&H EoS). Since 
only about one third of the ejected matter has a value very 
close to 0.5, the mass of ejected ®®Ni is certainly smaller than 
~0.015Mo. The values in our models are higher than those 
in previous simulations of SNe from O-Ne-Mg cores (Mayle & 
Wilson 1988). This points to important differences in the neu¬ 
trino treatment. Spectral Boltzmann transport calculations have 
recently found early ejecta with Y^ around 0.5 and higher also 
in (artificial) explosions of more massive progenitors (Buras 
et al. 2005, Erohlich et al. 2005). The reason for this differ¬ 
ence compared to the older models is a refined description of 
neutrino spectra formation and in particular of charged-current 
neutrino-nucleon interactions, including the weak magnetism 
corrections that were pointed out to be relevant by Horowitz 
( 2002 ). 

4. Conclusions 

Our ID simulations of SN explosions from the collapse of O- 
Ne-Mg cores suggest that such SNe are powered by neutrino 
heating and by the neutrino-driven wind of the newly formed 
NS, similar to what Woosley & Baron (1992) found in case of 
the accretion-induced collapse of white dwarfs to NSs. Such 
events have a low explosion energy (~ 0.1x10®' erg) and pro¬ 
duce little ®®Ni (:$ 10“^ Mq). Most of the ejecta expand initially 
with velocities of 2^xl0"'km s“', a small fraction has nearly 
10®kms *. This is significantly faster than in SNe of more 
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massive progenitors. Of course, sweeping up the matter of the 
stellar mantle and envelope, the shock will decelerate, and the 
ejecta velocities after shock breakout from the stellar surface 
will be correspondingly lower. During the first second of the 
explosion, the ejected matter has 0.46 % 0.53 and modest 

entropies (10 % sKk^jhy) $ 40). Such conditions exclude that 
r-process elements are formed in this matter during this early 
phase of the explosion of O-Ne-Mg cores: is too large for 

the “classical” low-entropy r-process and s is too low for high- 
entropy r-processing. The ejecta in our models, however, do 
not show the vast overproduction of some very neutron-rich, 
rare isotopes like *^Kr, which was made in low-Tg material 
(Te ^ 0.44) in previous simulations (Mayle & Wilson 1988), 
and which was interpreted as a severe constraint to the rate 
of such events. Our models yield considerably less energetic 
explosions than previous simulations and show other signifi¬ 
cant differences in the dynamics and explosion characteristics. 
These are probably mainly linked to the improved treatment of 
neutrino transport and neutrino-matter interactions. 

The small explosion energy obtained in our simulations is 
more consistent than previous explosion models of O-Ne-Mg 
cores with the low present expansion velocities (~ 1500 km s“') 
of the filaments of the Crab remnant of SN 1054 (Davidson 
& Fesen 1985), corresponding to a low kinetic energy of 0.6- 
1.5x10^^’erg for an ejecta mass of 4.6 + 1.8 M© in ionized and 
neutral gas (Fesen et al. 1997). The energy could be even some¬ 
what larger if there were several solar masses of material in an 
undetected, extended halo. 

While our simulations are spherically symmetric, we do not 
expect any qualitative changes in the multi-dimensional case, 
and probably only a modest increase of the explosion energy. 
Since very fast outflow develops on a relatively short timescale 
after core bounce, nonradial hydrodynamic instabilities are un¬ 
likely to have time to merge and grow to very large structures or 
global asymmetry before the anisotropic pattern freezes out in 
the accelerating expansion (Scheck et al. 2006, in preparation). 
Therefore the recoil velocity of the NS due to anisotropic mass 
ejection should remain fairly small (see Scheck et al. 2004), 
in agreement with speculations by Podsiadlowski et al. (2004). 
Corresponding 2D simulations are in progress. 

Our models have also an important bearing on the nuclear 
EoS constraints deduced by Podsiadlowski et al. (2005) from 
the low-mass Pulsar B of the double pulsar J0737-3039, which 
has a gravitational mass of Me - 1.249 + 0.001 M©. Provided 
the progenitor model we use is valid, the mass loss of the col¬ 
lapsing O-Ne-Mg core during the explosion leaves the neutron 
star with a baryonic mass of Mq = 1.36 + 0.002 M©. The error 
range accounts approximately for variations associated with the 
employed EoS and the wind ablation after our simulations are 
terminated. Our value implies a systematic left shift and reduc¬ 
tion of the “acceptance rectangle” in Eig. 3 of Podsiadlowski et 
al. (2005). Combined with the recent measurement of a pulsar 
of 2.1 +0.2 M© in PSR J0751 + 1807 (Nice et al. 2005), which is 
the largest well determined NS mass so far, this lends viability 
only to a limited number of NS EoSs which allow for a suf¬ 
ficiently large maximum mass and whose Mg/M o-curves pass 
through the acceptance rectangle. 


Based on our findings one might speculate that (B-) stars 
around OM© are also the progenitors of some of the sublu- 
minous Type II-P supernovae mentioned in the introduction. 
In fact, their peculiarities would be explained in a very natu¬ 
ral way. The low peak luminosity and extended plateau phase 
could result from the combination of a low hydrogen envelope 
mass (^ 6 M©) with low expansion velocities (^ 3000 km/s). 
The small mass of radioactive ^®Ni would explain the low tail- 
luminosity of these objects. An alternative interpretation of 
subluminous Type II-P supernovae is the explosion of rather 
massive stars with extended envelopes, but otherwise more 
“normal” explosion energies. This connection is supported by 
the long duration of the plateau phase of many of these events 
(cf. Pastorello et al. 2004, and refs, therein). Provided the ob¬ 
servations cover the full duration of the plateau, SN 1997D 
and 2003gd may still be viable cases for explosions of stars 
with main sequence masses around 9 M© (Hendry et al. 2005). 
There is a clear difference between such stars and more massive 
supernova progenitors. The former eject very little amounts 
(some 10“^ M©?) of oxygen only, like SN 1997D (Chugai & 
Utrobin 2000), whereas the latter produce up to a solar mass or 
more. Therefore one might be able to distinguish between the 
two scenarios on the basis of observations by measuring the 
oxygen lines in the late nebular spectra. 

Acknowledgements. We are grateful to K. Nomoto for pro¬ 
viding us with the initial data and for encouraging discussions. 
We thank R. Buras and M. Rampp for their input to this work, 
A. Marek especially for the W&H nuclear EoS, and P. Mazzali, 
A. Pastorello, D. Sauer and V. Utrobin for very helpful dis¬ 
cussions about SN observations. This work was supported by 
the Sonderforschungsbereich 375 “Astroparticle Physics” of 
the Deutsche Eorschungsgemeinschaft and by the International 
Max Planck Research School (IMPRS). Supercomputer time at 
the Rechenzentrum Garching is acknowledged. 


References 

Arnett, W.D. 1975, ApJ, 195, 727 

Baron, E., Cooperstein, J., & Kahana, S. 1987, ApJ, 320, 300 
Bethe, H.A. &. Wilson, J.R. 1985, ApJ, 295, 14 
Buras, R., Rampp, M., Janka, H.-T, & Kifonidis, K. 2005, A&A, in 
the press, astro-ph/0507135 
Burrows, A. & Lattimer, J. 1985, ApJ, 299, L19 
Chugai, N.N. & Utrobin, V.P. 2000, A&A, 354, 557 
Davidson, K. & Fesen, R.A. 1985, ARAA, 23, 119 
Davidson, K. et al. 1982, ApJ, 253, 696 
Fesen, R.A., Shull, J.M., & Hull, A.P. 1997, AJ, 113, 354 
Frbhlich, C. et al. 2004, ApJ, in the press, astro-ph/0410208 
Fryer, C.L., Benz, W., Herant, M., & Colgate, S. 1999, ApJ, 516, 892 
Gott, J.R., Gunn, J.E., & Ostriker, J. P. 1970, ApJ, 160, L91 
Hamuy, M. 2001, ApJ, 582, 905 

Hartmann, D., Woosley, S.E., & El Eid, M.E 1985, ApJ, 297, 837 
Hashimoto, M., Iwamoto, K., & Nomoto, K. 1993, ApJ, 414, L105 
Hendry, M.A. et al. 2005, MNRAS, 359, 906 
Hillebrandt, W. 1978, Space Science Reviews, 21, 639 
Hillebrandt, W. 1982, A&A, 110, L3 

Hillebrandt, W., Nomoto, K., & Wolff, R.G. 1984, A&A, 133, 175 
Horowitz, C.J. 2002, Phys. Rev. D, 65, 043001 
Janka, H.-T. 1999, unpublished 



6 


RS. Kitaura et. al.: Explosions of 0-Ne-Mg cores and the Crab supernova 


Langanke, K., Martinez-Pinedo, G., Sampaio, J.M., et al. 2003, 
Phys. Rev. Lett., 90, 241102 
Lattimer, J. & Swesty, R 1991, Nucl. Phys. A, 535, 331 
Marek, A., Dimmelmeier, H., Janka, H.-T., Muller, E., & Buras, R. 

2005, A&A, in the press, astro-ph/0502161 
Mayle, R. & Wilson, J.R. 1988, ApJ, 334, 909 
Miyaji, S. & Nomoto, K. 1987, ApJ, 318, 307 
Miyaji, S., Nomoto, K., Yokoi, K., & Sugimoto, D. 1980, 
Astron.Soc.Jap., 32, 303 

Nice, D.J., Splaver, E.M., Stairs, I.H., et al. 2005, ApJ, submitted, 
astro-ph/0508050 

Nomoto, K. 1983, Int.Astron.U., 101, 139N 
Nomoto, K. 1984, ApJ, 277, 791 
Nomoto, K. 1987, ApJ, 322, 206 

Nomoto, K., Sparks, W.M., Eesen, R.A., et al. 1982, Nature, 299, 803 
Pastorello, A. et al. 2004, MNRAS, 347, 74 
Pastorello, A. et al. 2005, MNRAS, submitted 
Podsiadlowski, Ph., Langer, N., Poelarends, A.J.T., et al. 2004, ApJ, 
612, 1044 

Podsiadlowski, Ph., Dewi, J.D.M., Poelarends, R, et al. 2005, 
MNRAS, 361, 1243 

Rampp, M. & Janka, H.-T. 2002, A&A, 396, 361 
Scheck, L., Plewa, T, Janka, H.-T, et al. 2004, PRL, 92, 011103 
Sumiyoshi, K., Terasawa, M., Mathews, G.J., et al. 2001, ApJ, 562, 
880 

Takahara, M., Hino, T., Oda, T., et al. 1989, Nucl. Phys. A, 504, 167 
Wanajo, S., Tamamura, M., Naoki, 1., et al. 2003, ApJ, 593, 968 
Wheeler, J.C., Cowan, J.J., & Hillebrandt, W. 1997, ApJ, 493, LlOl 
Woosley, S.E. & Baron, E. 1992, ApJ, 391, 228 
Woosley, S.E., Weaver, T.A., & Taam, R. 1980, in Type I Supernovae, 
ed. J.C. Wheeler (Austin: University of Texas), p. 96 
Zampieri, L. et al. 2003, MNRAS, 338, 711 



